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Abstract 

We present a detailed study of the total pp cross section for scalar Higgs pro- 
duction to next-to-next-to leading order in a s at LHC energies, and of the pp cross 
section at the Tevatron, combining an implementation of the solutions of the par- 
ton evolution equations at the 3-loop order with the corresponding hard scatterings, 
evaluated at the same perturbative order. Our solutions of the DGLAP equations 
are implemented directly in x-space and allow the study of the dependence of the 
results on the factorization (/ip) and renormalization scales (/zr) typical of a given 
process, together with the stability of the perturbative expansion. The input sets for 
the parton evolutions are those given by Martin, Roberts, Stirling and Thorne and 
by Alekhin. Results for K-factors are also presented. The NNLO corrections can be 
quite sizeable at typical collider energies. The stability region of the perturbative 
expansion is found when fiR > mjj ~ uf- 
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Figure 1: The leading order diagram for Higgs production by gluon fusion 

1 Introduction 

The validity of the mechanism of mass generation in the Standard Model will be tested 
at the new collider, the LHC. For this we require precision studies in the Higgs sector 
to confirm its existence. This program involves a rather complex analysis of the QCD 
backgrounds with the corresponding radiative corrections fully taken into account. Studies 
of these corrections for specific processes have been performed by various groups, to an 
accuracy which has reached the next-to-next-to-leading order (NNLO) level in a s , the 
QCD coupling constant. The quantification of the impact of these corrections requires 
the determination of the hard scattering partonic cross sections up to order af, together 
with the DGLAP kernels controlling the evolution of the parton distributions determined 
at the same perturbative order. Therefore, the study of the evolution of the parton 
distributions, using the three-loop results on the anomalous dimensions [T], is critical for 
the success of this program. Originally NNLO predictions for some particular processes 
such as total cross sections [2j have been obtained using the approximate expressions for 
these kernels [3J. The completion of the exact computation of the NNLO DGLAP kernels 
motivates more detailed studies of the same observables based on these exact kernels and 
the investigation of the factorization and renormalization scale dependences of the result, 
which are still missing. In this work we are going to reanalyse these issues from a broader 
perspective. Our analysis is here exemplified in the case of the total cross sections at 
the LHC (pp) and at the Tevatron (pp) for Higgs production using the hard scatterings 
computed in [2] and their dependence on the factorization and renormalization scales. The 
DGLAP equation is solved directly in x-space using a method which is briefly illustrated 
below and which is accurate up to order a 2 s . Our input distributions at a small scale 
will be specified below. We also analyse the corresponding K-factors and the region of 
stability of the perturbative expansion by studying their variation under changes in all 
the relevant scales. It is shown that the NNLO corrections are sizeable while the region 
of reduced scale dependence is near the value tuh = with [Mr around the same value 
but slightly higher. 

2 Higgs production at LHC 

The Higgs field, being responsible for the mechanism of mass generation, can be radiated 
off by any massive state and its coupling is proportional to the mass of the same state. 
At the LHC one of the golden plated modes to search for the Higgs is its production 
via the mechanism of gluon fusion. The leading order contribution is shown in Fig. QJ) 
which shows that dependence of the amplitude is through the quark loop. Most of the 
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Figure 2: A typical NLO diagram for Higgs production by gluon fusion 



contribution comes from the top quark, since this is the heaviest quark and has the largest 
coupling to the Higgs field. NLO and NNLO corrections have been computed in the last 
few years by various groups. A typical NLO correction is shown in Fig. (J2J). In the infinite 
mass limit of the quark mass in the loop (see |H for a review), an effective description of 
the process is obtained in leading order by the Lagrangian density 
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being the contribution of N H heavy fermion loops to the QCD beta function. This effective 
Lagrangian can be used to compute the radiative corrections in the gluon sector. A 
discussion of the NNLO approach to the computation of the gluon fusion contributions to 
Higgs production has been presented in [2J, to which we refer for more details. We recall 
that in this paper the authors presented a study for both scalar and pseudoscalar Higgs 
production, the pseudoscalar appearing in 2-Higgs doublets models. The diagonalization 
of the mass matrix for the Higgs at the minimum introduces scalar and pseudoscalar 
interactions between the various Higgs and the quarks, as shown from the structure of 
the operator 2 below in eq. JHJ). In the large top-quark mass limit the Feynman rules 
for scalar Higgs production (H) can be derived from the effective Lagrangian density [H] 



Cf ff = G R $ H (x) O (x) with 0{x) = - - Gl v (x) G a ^ v (x) , 
whereas the production of a pseudo-scalar Higgs (A) is obtained from 
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where $ H (x) and $ A (x) represent the scalar and pseudo-scalar fields respectively and nf 
denotes the number of light flavours. G£ u is the field strength of QCD and the quark 
fields are denoted by g^. We refer the reader to [2] for further details. 

Using the effective Lagrangian one can calculate the total cross section of the reaction 



H 1 (P 1 )+H 2 (P 2 )^E + X, (4) 

where Hi and H 2 denote the incoming hadrons and X represents an inclusive hadronic 
state and B denotes the scalar or the psudoscalar particle produced in the reaction. The 
total cross section is given by 

Otot = p /Arf B 7\ [ dx 1 [ dx 2 f a (Xl, [l 2 ) f b (x 2 , [J 2 ) 
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with x = — , S=(P 1 + P 2 ) 2 , (5) 

where the factor 1/(N 2 — 1) is due to the average over colour. The parton distributions 
f a {y,^ 2 ) (a,b = q,q,g) depend on the mass factorization/renormalization scale \i. A^b 
denotes the partonic hard scattering coefficient computed with NNLO accuracy. 



3 The NNLO Evolution 

We summarize the main features of the NNLO DGLAP evolution. As usual we introduce 
singlet (+) and non-singlet (— ) parton distributions 

n f 

q^ = qi ±q^ q ^=Y^q? ] (6) 
whose evolution is determined by the corresponding equations 



d ( q^(x,Q 2 ) \ _ ( P qq (x,a s m) P qg (x, a s (Q 2 )) \ / q^(^Q 2 ) \ m 
dlogQ^ g ( x ,Qi) ) \P gq (x,a s (Q 2 )) P gg (x,a s (Q 2 )) J { g(x,Q 2 ) J [ ' ] 

for the singlet combination and a scalar one for the non-singlet case 

-^A^q\-\x,Q 2 ) = P NS (x,a s (Q 2 ))®q^(x,Q 2 ). (8) 
The convolution product is defined by 

We recall that the perturbative expansion, up to NNLO, of the kernels is 

P(x, a s ) = a s P (0) (x) + a 2 s P {1) (x) + a 3 s P {2) (x) + . . . . (10) 
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where a s = a s /(4ir). In order to solve the evolution equations directly in x-space (see |H] 
for an NLO implementation of the method), we assume solutions of the form [Zj 

f( x n 2 ) = V An ^ log n + a (O 2 ) V Bn ^ log" Qs ^^ 

i 2/ / n2\ C n (x) CL S (Q ) M n 

+a ' W) £~' og S^I) 

for each parton distribution /, where Qo defines the initial evolution scale. The ansatz 
is introduced into the evolution equations and used to derive recurrence relations for its 
unknown coefficients A n , B n , C n , involving polylogarithmic functions (HIEI which are then 
implemented numerically. 

This ansatz corresponds to a solution of the DGLAP equation accurate up to order 
a 2 (truncated solution). It can be shown that this ansatz reproduces the solution of 
the DGLAP equation in (Mellin) moment space obtained with the same accuracy in a s . 
Modifications of this ansatz also allow to obtain the so-called "exact" solutions of the 
equations for the moments [10J. These second solutions include higher order terms in a s 
and can be identified only in the non-singlet case. Exact approaches also include an exact 
solution of the renormalization group equation for the /3-function, which embodies the 
effects of the coefficients /3 , Pi and /3 2 to higher order in a s . The term "exact" is, however, 
a misnomer since the accuracy of the solution is limited to the knowledge of the first three 
contributions to the expansion in the beta function and in the kernels. It can be shown 
both for exact and for the truncated solutions that solving the equations by an ansatz in 
x-space is completely equivalent to searching for the solution in moment space, since in 
moment space the recursion relations can be solved exactly 0. 



4 Renormalization scale dependence 

For a better determination of the dependence of the perturbative cross section on the 
scales of a certain process it is important to keep these scales independent and study the 
behaviour of the corresponding hadronic cross section under their variation. In our case 
the two relevant scales are the factorization scale /if and the renormalization scale /i^ 
which can be both included in the evolution by a rearrangement of the evolution kernels 
up to NNLO. 

The study of the dependence of the solution upon the various scales is then performed 
in great generality and includes also the logarithmic contributions log(/ii?//i/j) coming 
from the hard scatterings given in [2], where, however, only the specific point /ip = Hr = 
m H was considered. The separation of the scales should then appear not only in the hard 
scatterings but also in the evolution equations. This issue has been addressed in [TH] and 
can be reconsidered also from x-space [Zj using the x-space logarithmic ansatz (fH| . 

The scale dependence of the parton distribution functions is then expressed by a 
generalized DGLAP equation 

d 

d l nfJ 2 /4) = Pijfa vl, lAi) ® fj( x > /4> /4) > ( 12 ) 

where (ip is now a generic factorization scale. 

Generally speaking, both the kernels and the PDF's have a dependence on the scales 
Hf and /iR, and formally, a comparison between these scale is always possible up to a 
fixed order by using the renormalization group equations for the running coupling a s . 
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The renormalization scale dependence of the ansatz (jTT| that solves (fT2~| is obtained 
quite straighforwardly by a Taylor expansion of the running coupling a s (fMp) in terms of 
a s (/4) 



a 9 (fip) = a s ([f R ) 



4tt (4tt) 2 



-(3 2 L 2 + faL) 



(13) 



where the dependence is included in the factor L = ]n(jip / /jl r ) , and the coefficients of 
the /3-function, (the fa) are listed below 
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As usual we have set 

A c = 3, 
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T f = T R n f = -n f , 



(14) 



(15) 



where Nc is the number of colors and rif is the number of active flavors. This number is 
varied as we step into a region characterized by an evolution scale // larger than a specific 
quark mass (// > m q ). Also the NNLO matching conditions across flavor thresholds [TT] . 
|12j are implemented. 

Since the perturbative expansion of eq. (fTU|) contains powers of a s (fi F ) which can be 
related to the value of oc s (/i R ) by (fT3|) . from 
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substituting eq. (fT3|) into (jUJ), we obtain the corresponding expression of the kernels 
organized in powers of a s (ji R ) up to NNLO, and it reads (TOJ 
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The implementation of the method in x-space is quite straightforward and allows us 
to perform a separate study of the predictions in terms of fip and fiR. 



5 Numerical Results 

The use of the NNLO evolution of the parton distributions together with the results of [2] 
allows us to provide accurate predictions for the total cross section for Higgs production. 
Here we summarize and discuss our numerical results. 
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We use as inital conditions at low scales the sets of distributions given by MRST (THJ 
and Alekhin [Hj. Our final plots refer to center-of-mass energies which are reachable at 
the LHC, with 14 TeV being the largest one achievable in a not so distant future, and 
at the Tevatron, where we have selected the corresponding value as 2 TeV. We have also 
taken the Higgs mass run as a parameter in the prediction, with an interval of variability 
which goes from a light to a heavy Higgs (100 GeV to 300 GeV). Therefore //p, and mg 
are studied chosing various combinations of their possible values in the determination of 
total cross sections at leading (&lo), next-to-leading (<Jnlo), an d next-to-next-to-leading 
order (unnlo)- We present both standard two-dimensional plots and also some three 
dimensional plots in order to characterize in detail the structure of the region of stability 
of the perturbative expansion. We have also evaluated the K-factors for the total Higgs 
cross section at NLO, defined by 

K, = *™ (18) 

and at NNLO 

K 2 = (19) 

° NLO 

The study of the K-factors has been performed first by keeping the three scales equal 
(/xf = Hr = rnn) and then letting them vary around the typical value van. A second 
set of studies has been performed by taking typical values of ran and varying the value 
of the renormalization scale. We present in Figs. 01 LO and NLO results for the total 
Higgs cross sections at 2 TeV and at 14 TeV. The LO cross sections increase by a factor 
of approximately 100 as we change the energy from 2 TeV to 14 TeV (Figs. [3-a and 03-d) 
and sharply decrease as we raise the mass of the Higgs boson. At 14 TeV the range of 
variation of o^o is between 30 and 5 pb, with the highest value reached for m# = 100 
GeV. 

In the same figure we compare LO, NLO and NNLO cross sections at these two typical 
energies. It is quite evident that the role of the NLO corrections is to increase by a factor 
of approximately 2 the LO cross section bringing the interval of variation of Onlo between 
60 and 10 pb, for an increasing value of m#. NNLO corrections at 14 TeV increase these 
values by an additional 10 per cent compared to the NLO prediction, with a growth 
which is more pronounced for the set proposed by Alekhin. The two sets give coincident 
predictions for values of m H > 200 GeV. At lower values of m# (m H = 100 GeV) the 
NNLO predictions for a are quite different in the two sets. At the highest energy the 
cross sections evalutated in both models remain below 80 pb, while at Tevatron energies 
they are around 2 pb. 



5.1 K-factors 

A precise indication on the impact of the NLO/NNLO corrections and the stability of the 
perturbative expansion comes from a study of the K-factors K 1 and K 2 , defined above. 
From the plots in Figure SI the different behaviour of the predictions derived from the two 
models for the parton distributions is quite evident. At 14 TeV NLO K-factors from both 
models are large, as expected, since the LO dependence is not indicative of perturbation 
theory. The increase of ctnlo compared to the LO predictions is between 70 and 95 per 
cent. 

The two models also show a quite different behaviour for an increasing va#: the impact 
of the NLO corrections to the LO result predicted by Alekhin decrease for an increasing 
m H , while the MRST model has the opposite behaviour. The trend of the Alekhin set 
is reversed when we move to NNLO. The two K 2 factors grow between 1.05 and 1.35 at 
14 TeV, with the corrections predicted by Alekhin being the largest ones. The evaluation 
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of the overall impact of this growth on the size of these corrections should, however, also 
keep into consideration the fact that these corrections are enhanced in a region where the 
cross section is sharply decreasing (Figs.EJ. At 2 TeV K x is quite stable as m H increases 
(MRST set), varying between 2 and 2.2. In Alekhin's model this K-factor is smaller and 
its range is between 1.7 and 2. At the Tevatron the NNLO corrections are quite sizeable, 
with a factor K 2 between 1.4 and 1.8, as we increase the mass of the Higgs scalar. 

5.2 Renormalization/factorization scale dependence 

Now we turn to an analysis of the dependence of our results on /ip and fiR. In Figs. El 
we perform this study by computing a as a function of the Higgs mass for an incoming 
energy of 2 and 14 TeV and choose 

/4 = = 2m H . (20) 

We have seen that for a typical Higgs mass around 100 GeV and 14 TeV of energy (for 
m H — — Hr) the cross section doubles when we move from LO to NNLO, and a 
similar trend is also apparent if we fix the relations among the scales as in (f20|) . In this 
case, however, the impact of the NLO and NNLO corrections is smaller, a trend which 
is apparently uniform over the whole range of the Higgs mass explored. For m# = 100 
GeV the scalar cross section ctnnlo is around 58 pb for coincident scales, while a different 
choice, such as ([271 lowers it to approximately 45 pb. At Tevatron energies the variations 
of the cross section with the changes of the various scales are also sizeable. In this case 
for m H = 100 GeV the LO, NLO and NNLO predictions (0.6, 1.2 and 1.6 pb respectively) 
change approximately by 10-20 per cent if we include variations of the other scales as 
well. A parallel view of this trend comes from the study of the dependence of the K- 
factors. This study is presented in Figs. El The interval of variation of the K-factors 
is substantially the same as for coincident scales, though the trends of the two models 
[T8] and [14] is structurally quite different at NLO and at NNLO, with several cross-overs 
among the corresponding curves taking place for m# around 200 GeV. Another important 
point is that the values of K 2 are, of course smaller than K\ over all the regions explored, 
signaling an overall stability of the perturbative expansion. 

5.3 Stability and the Choice of the Scales 

The issue of determining the best of values of m# /ip and \ir in the prediction of the 
total cross section is a rather important one for Higgs searches at LHC. We have therefore 
detailed in Figs. and [HI a study of the behaviour of our results varying the renormalization 
scale hr at a fixed value of the ratio between fip and m H . In these figures we have chosen 
two values for the ratio between these two scales. Apart from the LO behaviour of the 
scalar cross section, which is clearly strongly dependent on the variation of both scales (see 
Figs, (a) and (d)), and does not show any sign of stability since the cross section can be 
drastically lowered by a different choice of fip, both the NLO and the NNLO predictions 
show instead a clear region of local stability for fiR > \ip but not too far away from the 
"coincidence region" \ip = fiR = m H . This can be illustrated more simply using Fig. E^b) 
as an example, where we have set the incoming energy of the p-p collision at 14 TeV. In 
this case, for instance, we have chosen ran — Hf — 100 GeV {C = 1), and it is clear 
from the plots that a plateau is present in the region of hr ~ 130 GeV. Similar trends are 
also clearly visible at NNLO, though the region of the plateau for the scalar cross section 
is slightly wider. Also in this case it is found that the condition hr > fip generates a 
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reduced scale dependence. Away from this region the predictions show a systematic scale 
dependence, as shown also for the choice of C — 1/2 in the remaining figures. In Figs. [HI 
we repeat the same study for the K-factors, relaxing the condition on the coincidence 
of all the scales and plotting the variations of K 1 and K 2 in terms of In the case 
nT>H — — 100 GeV the plateau is reached for hr ~ 150 GeV for K\ and hr ~ 200 
GeV for K 2 . In the first case the NLO corrections amount to an increase by 100 per cent 
compared to the LO result, while the NNLO corrections modify the NLO estimates by 
about 20 per cent (MRST). Similar results are obtained also for fip = 50 GeV. In this 
case, at the plateau, the NLO corrections are still approximately 100 per cent compared 
to the LO result and the NNLO corrections increase this value by around 15 per cent 
(MRST). 

5.4 Energy Dependence 

The energy dependence of the NNLO predictions for the total cross sections and the 
corresponding K-factors at the LHC are shown in Figs. ITllfTBl where we have varied the 
ratio C = fiF/rrin and k = jj? R / ji 2 F in order to illustrate the variation of the results. 
The cross sections raises sharply with energy and the impact of the NNLO corrections 
is significant. The K-factors, in most of the configurations chosen, vary between 1 and 
2.2. We have chosen the MRST input. The behaviour of the K-factors is influenced 
significantly by the choice of the ratio (k) between /i# and ^f- In particular, in Figs. fT5l 
the NNLO K-factors increase with \fS for k = 2 , the center of mass energy, which is not 
found for other choices of scales. The case k — 1/2 is close in behaviour to the coincident 
case n R = fip. The overall stability of the K-factors is clearly obtained with the choice 
k — 1. We have finally included in a set of tables our numerical predictions in order to 
make them available to the experimental collaborations. 

6 Conclusions 

A study of the NNLO corrections to the cross section for Higgs production has been 
presented. We have inplemented the exact three-loop splitting functions in our own 
parton evolution code. We used as initial conditions (at small scales) the boundary values 
of Martin, Roberts, Thorne and Stirling and of Alekhin. This study shows that the impact 
of these corrections are important for the discovery of the Higgs and for a reconstruction of 
its mass. The condition of stability of the perturbative expansion is also quite evident from 
these studies and suggests that the optimal choice to fix the arbitrary scales of the theory 
are near the coincidence point, with /i# in the region of a plateau. The determination of 
the plateau has been performed by introducing in the perturbative expansion and in the 
evolution a new independent scale (/zr) , whose variation allows to accurately characterize 
the properties of the expansion in a direct way. 

While this paper was being completed several authors have presented studies of the 
total Higgs cross section based on threshold resummation of soft or soft-plus-virtual log- 
arithms, see (THJ , [El, |T7j . Our work is based on exact NNLO partonic cross sections. 
Another relevant paper which recently appeared is |18j . 
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Figure 3: Cross sections for scalar Higgs production at the LHC and at the Tevatron with 
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Figure 4: K-factors for scalar Higgs at NNLO/NLO and NLO/LO with \i\ = /4 
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Figure 6: K-factors for scalar Higgs at NNLO/NLO and NLO/LO with ii\ = (l/2)/4 
and ix p = 2m H 
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Figure 8: K-factors for scalar Higgs production at LHC, NNLO/NLO and NLO/LO as a 
function of /zr, with jip = Cm H and tuh = 100 GeV 
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Figure 9: Three-dimensional graphs for the LO cross sections for scalar Higgs production 
at LHC, as a function of [Ir and with /xp with a fixed value of m# 
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Figure 10: Three-dimensional graphs for the NLO cross sections for scalar Higgs produc- 
tion at LHC, as a function of hr and with \ip with a fixed value of tuh 
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Figure 11: Three-dimensional graphs for the NNLO cross sections for scalar Higgs pro- 
duction at LHC, as a function of /ir and with \ip with a fixed value of mu 



19 




(a) m H = HOGeV 

Figure 12: Three-dimensional graphs for the ^-factor (Jnlo/^lo f° r scalar Higgs produc- 
tion at LHC, as a function of jiR and \ip and at a fixed value of m#. 
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Figure 13: Three-dimensional graphs for the i^-factor (Jnnlo/ '&nlo f° r scalar Higgs pro- 
duction at LHC, as a function of jiR and \ip and at a fixed value of m#. 
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Figure 14: Cross sections and ^-factors for scalar Higgs production at LHC as a function 
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Figure 15: Cross sections and K- factors for the scalar Higgs production at LHC as a 
function of with /ip — Cm H , with fip = k[i 2 R and mu = 114 GeV. MRST inputs have 
been used. 
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Figure 16: Cross sections and ^-factors for the scalar Higgs production at LHC as a 
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Table 1: Values of the cross sections and K- factors for scalar Higgs production at LHC 
as a function of \/S with \ip = mn, with \j? F = n 2 R and ran = 114 GeV. MRST inputs 
have been used. 
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Table 2: Values of the cross sections and K- factors for scalar Higgs production at LHC 
as a function of \fS with \ip = 2m#, with \j? F = n 2 R and ran = 114 GeV. MRST inputs 
have been used. 
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Table 3: Values of the cross sections and K- factors for scalar Higgs production at LHC as 
a function of \fS with \xf = (l/2)m#, with /ip = n 2 R and ran = 114 GeV. MRST inputs 
have been used. 
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Table 4: Values of the cross sections and ^-factors for scalar Higgs production at LHC 
as a function of \fS with fip — m H, with /i% = 2(i 2 R and mn = 114 GeV. MRST inputs 
have been used. 
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Table 5: Values of the cross sections and K- factors for scalar Higgs production at LHC 
as a function of \fS with /ip = ^ttih, with fj, F = 2/i R and m# = 114 GeV. MRST inputs 
have been used. 
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Table 6: Values of the cross sections and K- factors for scalar Higgs production at LHC as 
a function of v^S* with \ip = (l/2)m#, with pL% = 2\j? R and ran = 114 GeV. MRST inputs 
have been used. 
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Table 7: Values of the cross sections and ^-factors for scalar Higgs production at LHC as 
a function of \J~S with \xf = rriH, with pip = (l/2)/x|, and mn = 114 GeV. MRST inputs 
have been used. 
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Table 8: Values of the cross sections and K- factors for scalar Higgs production at LHC as 
a function of v^S* with \ip = 2m#, with fj,p = (1/2)//^ and ran = 114 GeV. MRST inputs 
have been used. 
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Table 9: Values of the cross sections and -KT-factors for scalar Higgs production at LHC 
as a function of VS with /i F = (l/2)m H , with /i 2 F = (l/2)/4 and m H = 114 GeV. MRST 
inputs have been used. 
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